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We report the application of the density-matrix renormalization group method to a spatially 
anisotropic two-dimensional Hubbard model at half-filling. We find a deconfinement transition 
induced by the transverse hopping parameter ty from an insulator to a metal. Therefore, if ty is 
fixed in the metallic phase, increasing the interaction U leads to a metal-to-insulator transition at 
a finite critical U. This is in contrast to the weak-coupling Hartree-Fock theory which predicts a 
nesting induced antiferromagnetic insulator for any U > Q. 
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The metal-insulator transition (MIT), also called the 
Mott transition [l| , is certainly one of the most difficult 
challenge facing condensed matter theorists. Hubbard 
[i^, in a pioneering work, introduced a simple one-band 
Hamiltonian which has only two parameters, t for the ki- 
netic energy of the electrons and U for the local electron- 
electron interactions. This model is at half-filling the 
model of reference for the MIT. In D = 1, Lieb and Wu 
Q obtained an exact solution by using the Bethe ansatz. 
The ground state is an insulator for any U/W > 0, where 
W is the band width. Thus, the MIT occurs at the criti- 
cal value {U/W)c = 0. In infinite dimensions, the dynam- 
ical mean- field theory (DMFT) 0, ^ predicts a critical 
point at {U/W)c w 1. 

The discovery of layered materials, where the motion 
of electrons driving the low energy physics is mostly con- 
fined in the layers, has raised great interest into the two- 
dimensional (2D) Hubbard model. The physics at large 
U/W ^ 1 is now understood, the charge excitations are 
gapped, the spin excitation are described by the Heisen- 
berg Hamiltonian which has long-range order (LRO) at 
T = 0. But for U/W ^ 1, the physics is still unclear. 
Our current knowledge about the weak-coupling regime is 
mostly drawn from the Hartree-Fock approximation and 
from quantum Monte Carlo (QMC) simulations 0]- 
The QMC results agree qualitatively with the Hartree- 
Fock prediction that the ground state is a Slater insulator 
for any U/W > 0. However, in most recent studies such 
as in Ref . [7|] , even though considerable progress has been 
achieved in reaching larger systems, in the weak U regime 
where the eventual gap is small, reliable extrapolations of 
the QMC data remain difficult to achieve. It would thus 
be preferable to apply finite size scaling for data analysis 
instead of relying on extrapolations. 

More recently, extensions of the DMFT which include 
non-local fiuctuations, the dynamical cluster approxima- 
tion (DC A) d or the cellular DMFT have been 
applied to the 2D Hubbard model. The focus in these 
studies have mostly been to discuss the nature of the MIT 
within the paramagnetic solution of the DMFT equa- 
tions. A systematic comparison of the possible ordered 
or disordered ground states as function of the cluster sizes 
is still lacking. Therefore, the issue as to whether or not 
quantum fluctuations destroy the Hartree-Fock solution 



in the half-filled 2D Hubbard model in the small U regime 
remains open. 

In this letter, we show that insight into this prob- 
lem can be gained by studying the quasi- ID Hubbard 
model. We apply the two-step density-matrix renormal- 
ization group (DMRG) [llj to an array of coupled Hub- 
bard chains. We find that there is a deconfinement tran- 
sition from the ID insulator toward a metallic phase as 
the transverse hopping ty is increased from 0. Hence, 
for a fixed ty in the metallic phase, there is a quantum 
phase transition (QPT) at a finite quantum critical point 
{QCP){U/W)c- This suggests, as seen in the limit of in- 
finite dimensions [4], that by freezing out the local time 
dynamics, the Hartree-Fock approximation is unable to 
account for the physics of the Hubbard model even for 
weak interactions. We note that Biermann and cowork- 
ers applied the chain-DMFT to the quasi-lD Hub- 
bard model. They also found a deconfinement transition. 
However, they were restricted to paramagnetic solutions. 
Hence, unlike our study, they could not tell whether or 
not their metallic solution is the true ground state of the 
quasi- ID Hubbard model. 

Let us briefiy describe the two-step DMRG method in- 
troduced for coupled chains in Ref.pjJ. This method is 
generalized here for coupled multi-leg ladders. For a sys- 
tem of spins or electrons on an anisotropic square lattice 
of dimensions L,j: x Ly, the Hamiltonian may be written 

as: H = H.ntra + Hinter ^ J2i f^i + gHinter, whcrc hi are 
the Hamiltonians of ID systems; <? ^ 1 is the transverse 
coupling. Hinter is of the same magnitude as Hintra- hi 
can represent a single chain, a two-leg ladder or even a 
multi-leg ladder as illustrated in Fig[T]for a spin system. 
The DMRG analysis of H is done in two steps. In the first 
step, the m2 lowest eigenfunctions and eigenvalues En 
from different charge-spin sectors of the ID Hamiltonian 
hi are obtained by applying the conventional ID DMRG 
algorithm [l3|. During this step, mi states are kept such 
that the lowest m2 states are accurately computed. In 
the second step, H is projected onto the tensor product 
of the (/)„'s. Since the resulting effective 2D Hamiltonian 
is now ID in the transverse direction, it can be studied 
using the conventional DMRG. Let for a given size, Eq 
and Eni2 be the lowest and the highest DMRG energies 
kept for hi. If 5 ^ AE — Em^ — Eq (in practice we 
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FIG. 1: Clusters used as building blocks for the 2D square 
lattice: (a) chains; (b) two-leg ladders; (c) four-leg ladders: 
Jx = 1 is the coupling along the legs; Jy is the inter-cluster 
coupling, and Jr is the coupling between the rungs for the 
ladders. 



set < 10), the two-step DMRG retains high accu- 
racy. But as Lx increases, AE decreases, for a fixed g, it 
would be impossible to retain accuracy for arbitrary large 
Lx. Thus, successfully performing finite size analysis for 
the g induced QPT will depend on the magnitude of the 
critical value of g and on the density of the ID spectrum. 

In order to avoid uncertainties related to extrapola- 
tions in the regime where the gap is very small, it is 
preferable to perform finite size scaling analysis. In the 
vicinity of the QCP the gap, A{g) oc {gc — gY , where v 
is the correlation length exponent ( Lorentz invariance is 
assumed), and gc is the value of the g at the QCP. The 
product LA, where L is the linear size of the system is 
given by a universal function, LA = f{C{gc — g)L^^^) 
14| . where C is independent of L. At the critical point, 
LA = /(O) is independent of L. All the curves A{g) 
for different sizes should converge aX g = gc- In practice 
however there are finite size effects. It is necessary to 
extrapolate the different crossing points in order to pre- 
cisely locate gc- Before studying the MIT in the Hubbard 
model, we first apply the two-step technique to the QPT 
in quantum antiferromagnets (AFM). There are two mo- 
tivations for this prelude on spin systems. First, well 
controlled QMC are available, this will allow to gauge 
the performance of the two-step DMRG. Second, the en- 
ergy scales of the different systems roughly span those of 
the Hubbard model when the interaction is varied. 

In coupled chains with S = 1, the ground state is 
known both in ID and 2D. In ID, the s yst ems has a spin 
gap, the Haldane gap As^i = 0.41049 [13]. The correla- 
tion length Cg=i ~ 6 [18|. In 2D, we know from the rig- 
orous result that the ground state has an LRO, 
thus it is gapless and ^5=1 — 00. Hence, there should 
be a QPT at some critical Jy from a disordered to and 
ordered ground state. QMC studies predict Jy « 0.04 



21| and J^ = 0.043648(8) ;22]. In the two-step study. 



we applied periodic boundary conditions in the x direc- 
tion and open boundary conditions in the y direction. 
We show in Figl2fa) the finite size behavior of the spin 
gap As. It shows that in agreement with QMC, Jy is 
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FIG. 2: As X Lx as function of Jy (a) for coupled 5* = 1 
chains, (b) for coupled S — 2 chains; (c) As as function of 
Lx at two characteristics values below (Jy = 0.02) and above 
{Jy = 0.06) the QCP for S = 1 chains; (d) 5 = i two-leg 
ladders at Jr = 0.4; (e) and four-leg ladders with Jr = 1; (f) 
Jy as function of L^^ in coupled four- leg ladders. 



accurate DMRG estimate of Jy can be obtained by lo- 
cating the crossing point of consecutive x Ly systems 
and extrapolating. But for S" = 1, an extrapolation was 
not necessary, for x Ly > 16 x 17, the data con- 
verged around Jy = 0.043613. We performed the same 
analysis for coupled S = 2 chains. For a single chain 
As=2 = 0.0876 fli and ^5=2 « 49 [3. The scaled gaps 
are displayed in FigHfb), we find J^ w 0.0007. It is in- 
teresting to note that the values of Jy for S = 1 and for 
S = 2 are roughly consistent with the Schwinger bosons 
prediction, Jy x « J^ [23|, where is the ID cor- 

2 = 1.5701 for = 1, 
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relation length. We find Jy x 
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that Jy X 



^ — 1.5707 for S = 2- These results suggest 
= ^Jx- In Figl2Jc) we show Ag for two 
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typical values in a 5 = 1 system, 
Jy — 0.06 > Jy. The extrapolated values are consistent 
with respectively gapped and gapless phases. 

In S* = i two-leg ladders, if the coupling between the 



located between Jy = 0.040 and Jy — 0.045. A more rungs is Jr = Jx = 1, QMC studies [22| predict that Jy 
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0.3. We could not study the QPT because, the condition 
Jy <^ AE is not fulfilled. If we reduce Jr enough, the 
two-step DMRG becomes applicable. This is achieved 
when Jj. < 0.5. For instance in Figl^fd), we display the 
finite size behavior of As for Jr — 0.4 and x Ly ranging 
from 10 X 12 to 30 x 32. For these systems, we find 



L xL=16xl7 



Jy = 0.0993. An alternative to reducing Jr for S ~ 

1 



systems in order to study the QPT in = ^ systems is 
increasing the number of legs. When the number of legs 
increases. As on the ladder decreases, the system is thus 
closer to criticality, therefore a smaller Jy can induce a 
QPT. For the four-leg ladder with Jr = 1, the QMC 
predicts Jy ~ 0.07 [231 . In Figmje), we show systems of 
four-leg ladders ranging from 12 x 12 to 44 x 44, it can 
be seen that Jy « 0.08. In Figl2][f), we plot the crossing 
points of 12 X 12 and 20 x 20 to 36 x 36 and 44 x 44 
systems respectively. This yields a better estimate of the 
QCP, J^ = 0.0742. 

Let us now consider the anisotropic Hubbard model 
defined by hopping parameters = ty and 
a local interaction U. The non-interacting single par- 
ticle energies are, Ck — —txcosik^x) — tyCos{kyy) where 
k = (kxjky). Since the Fermi surface is nested at the 
momentum qn, e(k -I- qn) = — e(k). In the Hartree- 
Fock approximation, the consequence of nesting is that 
the metallic state becomes unstable against the forma- 
tion of a gap A^^ m exp(^^) and spin-density wave 
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for all U and the charge gap 
exp(-^) for U/W < 1, and « U for 



LRO. 
Ai^ 

U/W ^ 1. There is no LRO in ID, the gap opening 
cannot be explained in the Hartree-Fock approximation. 
But it is generally believed that the Hartree-Fock approx- 
imation is at least qualitatively correct when ty ^ 0. 

However, there is a regime of the quasi-lD model 
where a simple physical argument shows the failure of 
the Hartree-Fock approximation. If ty ^ Aj-°, inter- 
chain motion is prohibited, the electrons are confined into 
the chains. For any U > 0, the system would remain a 
Mott insulator. There will be an LRO of the Heisenberg 
type because the spin degrees of freedom are gapless in 
ID, and the small ty would yield an effective exchange, 
Jy = t^l A\^ . This regime cannot be described by the 
simple Hartree-Fock theory. This shows that a strong 
coupling like behavior extends even for small U in the 
confined regime. This somewhat overlooked regime of 
the Hubbard model was discussed for two-coupled Hub- 
bard chains It is shown in Ref.|24| that if ty <€. A\^ 



the systems is equivalent to the Heisenberg two-leg spin 
ladder. This is also implicit in the chain dynamical mean- 
field theory study which predicted that at half-filling, a 
finite ty was necessary to deconfine the electrons in the 
transverse direction. Hence, if we increase ty from the 
ID Mott insulator ai ty ~Q there are three possibilities: 
(i) the system remains a Mott insulator; (ii) there is a 
crossover from a Mott insulator towards a Slater insula- 
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FIG. 3: (a) Ac x L^, as function of ty ak, U — 1; (b) the 
deconfinement transition point ty aXU = 1 extrapolated from 
crossing points of two consecutive x Ly; (c) Ac as function 
of l/La; at U = 1 for two characteristic values of ty, below 
{ty = 0.002) and above {ty = 0.05) ty-, Ac x as function of 
U in ID (d) and 2D (e); (f) Phase diagram of the quasi-lD 
Hubbard model in the {U,ty) plane. 



tor; (iii) there is a QPT towards a metallic phase. 

The two-step DMRG results shown below are consis- 
tent with the case (iii). In this study, we kept up to 
TOi = 384 and m2 = 96 respectively during the first 
and second steps; the maximum truncation error was less 
than 1 X 10~^ in the first step and 1 x 10"'' in the second 
step for all the parameters investigated. The bulk of our 
calculations was on lattice sizes x Ly ranging from 
12 X 13 to 32 X 33, J7 = 0, 0.5, 1, 1.5, 2, 2.5, 3, and 
at least 10 different ty chosen from ty ~ to ty Ri A^^ . 
In a few cases, 40 x 41 and 48 x 49 systems were also 
studied. In FigOl^a) we show the scaled Ac as function 
of ty for J7 = 1 . The scaled gap displays the typical be- 
havior seen for spin systems. The data for different sizes 
converge near ty = 0.007. The finite size behavior of ty 
(the crossing points of LxAc of consecutive systems) is 
shown in Fig[3l^b). The extrapolation yields ty = 0.0056 
for U = 1. In Fig[3]Jc), we display Ac as function of L^ at 
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two typical values of ty above and below . The extrap- 
olated gap is in agreement with the qualitative behavior 
of Lj-Ac- Hence displays a deconfinement transition 
from a ID Mott insulator to a 2D metallic phase. Using 
the chain-DMFT, Biermann and coworkers pjj] find that 
at U = 2.6, the charge correlation exponent Kp jumps 
from 0.02 at ty = 0.16 to 1.01 at = 0.28. This shows 
a deconfinement transition with 0.16 ^ ty ^ 0.28. We 
could not however see the ty induced QPT for U ~ 2.6. 
The relation ty ^ AE was satisfied only for very small 
systems. The maximum value for which we could study 
the ty induced QPT is [/ = 2, we find t^ = 0.0689. 
We obtain ty = 0.18 at [/ = 2. 6 by extrapolating from 



smaller values of U. Since in Ref.[l2] they were restricted 
to paramagnetic solutions, they could not rule out a pos- 
sible AFM ground state of the Slater type. 

Once in the deconfined regime, we can induce a MIT 
by increasing U. For this purpose, we set ty = 0.05. For 
this value of ty we know from the calculations above that 
for any U ^ 1.5, the system is in the metallic phase. 
We thus expect a MIT at some U between U = 1.5 and 
U = 2, because as seen above ty = 0.0689 at [/ = 2. 
In Figl3fd),(e), we show the scaled Ac as function of U 
respectively in ID and 2D. In the 2D case (FiglJl^e)), 
For U < 1.875, L^A^ decays when we increase as 
ioi U — until it reaches Uc ~ 1.8 where, it starts to 
increase. This is to be contrasted to the ID case shown in 
FiglSjd) where there is no regime where L^Ac decreases 
when is increased, which implies ?7c = as we know 
from the Lieb-Wu solution. 

In FiglHl^f), we show the phase diagram of the quasi-lD 
Hubbard model. The deconfinement transition occurs for 
small, [/ < 1, at tJJ « Aj^/4. In a recent experiment, 
Pashkin and coworkers 25] studied the infrared response 
of the quasi-lD Fabre sah (TA/TTFjaPFg. Though this 
compound is nominally three-quarter filled, the presence 
of small dimerization renders it effectively half-filled. At 
ambient pressure, {TAITTF)2PFq is a Mott insulator. 
When external pressure is applied, the interchain transfer 
integral grows exponentially and the Mott gap rapidly de- 
creases until a deconfinement transition is reached. The 
interchain transfer integral that induces the deconfine- 
ment transition is approximately half of the Mott gap. 
This result is in reasonable agreement with our predic- 
tion. 

To conclude, let us comment on the implication of our 
result on the isotropic case, ty = tx- For a given U , when 
ty > ty, the system enters the metallic phase. It should 
remain in the metallic phase up to ty — t^, because, as 
soon as ty > ty, there is no other obvious process that 
will drive the system to a another phase when ty is fur- 
ther increased. In Ref.fsj the DC A was applied to the 
isotropic 2D Hubbard model at half-filling. It was found. 



for cluster sizes up to Nc = 64, that down to U = 4 
the paramagnetic solutions remained gapped. Although 
no gap was found for [/ < 4, it was assumed that the 
gap would open for larger clusters which were not acces- 
sible. However, our result suggests that this assumption 
may not be true. We reexamined Ref.Q data, plotting 
Uc{Nc), where Uc{Nc) is the estimated finite cluster QCP, 
we find Uc ~ 3.1 when Nc — > oo. 
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